Some new static charged spheres 



S.D. Maharaj *, S. Thimkkanesh ^ 

Astrophysics and Cosmology Research Unit, School of Mathematical Sciences, 
University of KwaZulu-Natal, Private Bag X54001, Durban 4OOO, South Africa. 

o 
o 

^ ! Abstract 

, We present new exact solutions for the Einstein-Maxwell system in static spherically 

symmetric interior spacetimes. For a particular form of the gravitational potentials 
and the electric field intensity, it is possible to integrate the system in closed form. 
For specific parameter values it is possible to find new exact models for the Einstein- 
CJ I Maxwell system in terms of elementary functions. Our model includes a particular 

charged solution found previously; this suggests that our generalised solution could 
be used to describe a relativistic compact sphere. A physical analysis indicates that 
the solutions describe realistic matter distributions. 
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1 Introduction 



Exact solutions of tlie Einstein-Maxwell system of field equations, for spher- 
■ ically symmetric gravitational fields in static manifolds, are necessary to de- 

^ scribe charged compact spheres in relativistic astrophysics. The solutions to 

the field equations generated have a number of different applications in rel- 
ativistic stellar systems. It is for this reason that a number of investigations 
have been undertaken on the Einstein-Maxwell equations in recent times. A 
comprehensive review of exact solutions and criteria for physical admissability 
is provided by Ivanov [Tj. A general treatment of nonstatic spherically sym- 
metric solutions to the Einstein-Maxwell system, in the case of vanishing shear 
was, performed by Wafo Soh and Mahomed [2| using symmetry methods. The 
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uncharged case was considered by Wafo Soh and Mahomed [3] who show that 
all existing solutions arise because of the existence of a Noether point sym- 
metry; the physical relevance of the solutions was investigated by Feroze et al 
[1]. The matching of nonstatic charged perfect fluid spheres to the Reissner- 
Nordstrom exterior metric was pursued by Mahomed et al [5j who highlighted 
the role of the Bianchi identities in restricting the number of solutions. 

In this paper, we seek a new class of solutions to the Einstein-Maxwell sys- 
tem that satisfies the physical criteria. We attempt to perform a systematic 
series analysis to the coupled Einstein-Maxwell equations by choosing a ratio- 
nal form for one of the gravitational potentials and a particular form for the 
charged matter distribution. This approach produces a number of difference 
equations, which we demonstrate can be solved explicitly from first princi- 
ples. A similar approach was used by Thirukkanesh and Maharaj [6] , Maharaj 
and Thirukkanesh [7] and Komathiraj and Maharaj [8J. They obtained ex- 
act solutions by reducing the condition of pressure isotropy to a recurrence 
relation with real and rational coefficients which could be solved by math- 
ematical induction. In this way new mathematical and physical insights in 
the Einstein-Maxwell field equations were generated. An advantage of this 
approach is that we generate new solutions to the Einstein-Maxwell system 
which contain uncharged solutions found previously: Maharaj and Leach [5], 
Tikekar ^0], Durgapal and Bannerji [11] and John and Maharaj [12], amongst 
others. 

We first express the Einstein-Maxwell system of equations for static spherically 
symmetric line element as an equivalent system using the Durgapal and Ban- 
nerji [T^ transformation in Section 2. In Section 3, we choose specific forms 
for one of the gravitational potentials and the electric field intensity, which 
reduce the condition of pressure isotropy to a linear second order equation in 
the remaining gravitational potential. We integrate this generalised condition 
of isotropy equation using the method of Frobenius in Section 4. In general 
the solution will be given in terms of special functions. We demonstrate that 
it is possible to find two category of solutions in terms of elementary functions 
by placing certain restriction on the parameters. We regain known charged 
Einstein-Maxwell models from our general class of models in Section 5. In 
Section 6, we discuss the physical features of the solutions found, plot the 
matter variables, and show that our models are physically reasonable. 



2 The field equations 

We assume that the spacetime is spherically symmetric and static which is 
consistent with the study of charged compact objects in relativistic astro- 
physics. In Schwarzschild coordinates (t, r, 9, (p) the generic form of the line 
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element is given by 

ds' = -e^'^^Uf + e^^^'Ur^ + r\de^ + sin^ ed<P^). (1) 
The Einstein field equations for the line element ([1]) can be written as 
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for charged perfect fluids. The energy density p and the pressure p are mea- 
sured relative to the comoving fluid 4-velocity u"" = e~^5Q and primes denote 
differentiation with respect to the radial coordinate r. The quantities associ- 
ated with the electric field are E, the electric field intensity, and a, the proper 
charge density. In the system (l2al) - pd|) . we are using units where the coupling 
constant = 1 and the speed of light c = 1. This system of equations 
determines the behaviour of the gravitational field for a charged perfect fluid 
source. A different, but equivalent form of the field equations, can be found if 
we introduce the transformation 
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where the parameters A and C are arbitrary constants. Under the transfor- 
mation ([3]) the system (!2a|) - pdl) has the equivalent form 
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where dots denote differentiation with respect to x. This system of equations 
governs the behaviour of the gravitational field for a charged perfect fluid 
source. When = the Einstein-Maxwell equations (l4a!) - fl4dp reduce to the 
uncharged Einstein equations for a neutral fluid. Equation (l4c|) is called the 
generalised condition of pressure isotropy and is the fundamental equation 
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that needs to be integrated to demonstrate an exact solution to the Einstein- 
Maxwell system of equations (Ha!) - (14dp . 



3 Choosing potentials 



Our objective is to find a new class of solutions to the Einstein-Maxwell system 
by making explicit choices for the gravitational potential Z and the electric 
field intensity E. We make the choice for Z as 

where a and h are real constants. Note that the choice ([5]) ensures that the 
gravitational potential e^^ is regular and well behaved in the stellar interior 
for a wide range of values of the parameters a and h. In addition, when a; = 
then Z = \ which ensures that there is no singularity at the stellar centre. A 
special case of ([5]) was studied by Komathiraj and Maharaj [13]. The choice 
dS]) does produces charged and uncharged solutions which are necessary for 
constructing realistic stellar models. On substituting ^ in (Hcl) we obtain 



4(1 + axf{\ + hx)y + 2(1 + aa;)[6(l + ax) - 2(6 - a)\y 



+ 



by 



y = 0, (6) 



Cx 

which is a second order differential equation. 

The differential equation may be solved if a particular choice of the electric 
field intensity E is made. For our purpose we set 

E"^ aa{b — a)x , , 

~C " {l + hxf ' ^ ' 

where a is a constant. The electric field intensity specified in ([7]) vanishes at 
the centre of the sphere; it is continuous and bounded in the stellar interior for 
wide range of values of x. The quantity E"^ has positive values in the interior of 
star for relevant choices of the constants a, a and h. Therefore the form given 
in (171) is physically reasonable to study the behaviour of charged spheres. With 
the choice ([7]) we can express (E]) in the form 

4(1 + axf [6(1 + ax) - {b - a)]y + 2a(l + ax) [6(1 + ax) - 2(6 - a)] y 
+a{h — a){h — a — aa)y = 0. (8) 
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In ([8]) we assume that a 7^ and a 7^ 6 so that the electric field intensity is 
present. When a = there is no charge. 



4 Solutions 



To find the solution of the Einstein-Maxwell system we need to integrate 
the master equation ([8]). We consider two cases on the integration process: 
a = ^ - 1 and a ^ - 1. 

n. I n. 



A.l The case a = - — 1 

^ a 



In this case the differential equation ([8]) becomes 

2(1 + ax) [6(1 + ax) - {h - a)]y + a [6(1 + ax) - 2(6 - a)\y = 0. 
Equation ([9]) is easily integrable and the solution can be written as 
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where ci and C2 are constants of integration. Therefore, the solution of the 
Einstein-Maxwell system f l4al) - (l4dl) becomes 
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Observe that because of the restrictions a 



I and h ^ a the charged 



solution (|TT1) does not have an uncharged limit. Therefore this solution models 
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a sphere that is always charged and cannot attain a neutral state. Note that 
the solution (fTTj) is expressed in a simple form in terms of elementary functions 
which facilitates a physical analysis of the matter and gravitational variables. 



4.2 The case a ^ - — 1 

^ 'a 



With a 7^ ^ — 1, equation is difficult to solve. Consequently we introduce 
the transformation 

y={l + axyU{l + ax), (12) 

where f/ is a function of (1 + ax) and d is constant. With the help of f|T2l) . the 
differential equation ([8]) can be written as 

4(1 + ax f [6(1 + ax) - {b - a)] U 

+2(1 + ax) [b{4:d + 1)(1 + ax) - 2{2d + 1)(6 - a)] U 



+ 



2bd{2d - 1)(1 + ax) - (6 - a) - - 1 - a - 



U = 0. 



Note that there is substantial simplification if we take 

- - 1 - a = 4c/l 
a 

Then f[T^ becomes 
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where 6 7^ 0. We observe that the point 1 + ax = ^ is a regular singular 
point of the differential equation f|T^ . Therefore, the solution of the differential 
equation iHM can be written in the form of an infinite series by the method 
of Frobenius: 
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where Cj are the coefficients of the series and r is the constant. To complete the 
solution we need to find the coefficients q as well as the parameter r explicitly. 

The indicial equation determines the value of r from 

cor{2r - 3) = 0. 
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As Co 7^ we must have r = 0orr = 3/2. We can express the structure for 
the general coefficient Cj in terms of the leading coefficient Cq as 



[(p + r - 1) (2p + 2r + 4c^ - 3) + d{2d - 1)] 
li {p + r)i2p + 2r-3) ^'^^ 



where the conventional symbol H denotes multiplication. We can verify the 
result (fT6l) using mathematical induction. We can now generate two linearly 
independent solutions to (JT^ with the help of (ITSl) and ( fT6l) . For the parameter 
value r = 0, we obtain the ffist solution 
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For the parameter value r = 3/2, we obtain the second solution 
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Since the functions Ui and U2 are linearly independent we have found the 
general solution to (fT4|) . Therefore, the solutions to the differential equation 
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Thus the general solution to the differential equation ([6]) , for the choice of the 
electric field ([7]), is given by 



y{x) = Aiyi{x) + A2y2{x), (19) 



where Ai and A2 are arbitrary constants and d'^ = (^^ — 1 — aj / 4. From (1191) 
and fHa|) - fl4d| ). the exact solution of the Einstein-Maxwell system becomes 
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We believe that this is a new solution to the Einstein-Maxwell system. In gen- 
eral the models in ( l20l) cannot be expressed in terms of elementary functions 
as the series in ( |T71) and ( |T8|) do not terminate. Consequently the solution 
will be given in terms of special functions. Terminating series are possible for 
particular values of a and b as we show in the next section. 



5 Elementary functions 



It is possible to generate exact solutions in terms of elementary functions from 
the series in flTI?]) . This is possible for specific values of the parameters a, b and 
a so that the series (fT7j) and (fT8|) terminate. Consequently two categories of 
solutions are obtainable in terms of elementary functions by placing restric- 
tions on the quantity ^ — 1 — a. We can express the first category of solution, 
in terms of the variable r, as 
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where ^ — 1 — a = 4n^ relates the constants a, 6, a and n. The second category 
of sohition, in terms of the variable r, is given by 
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where - — 1 — a = 4:n{n — 1) + 1 relates the constants a, h, a and n. Thus we 
have extracted two classes of solutions (^T^ and fl221) to the Einstein-Maxwell 
system in terms of elementary functions from the infinite series solution fll9p . 
This class of solution can be expressed as combinations of polynomials and 
algebraic functions. The simple form of fl2T|) and fl22|) helps in the study of the 
physical features of the model. 

From our general classes of solutions fl2Tl) and fl22|) . it is possible to generate 
particular solutions found for charged stars previously . If we take 6=1 and 
K = then it is easy to verify that the equation ( 12T1) becomes 
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^ and D2 = (i_^^'-3/2 • Also, equation ([22]) 
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becomes 
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Thus we have 



2 ~ (l-a)"-2 • 

regained the Komathiraj and Maharaj [13] charged model; our solutions allow 
for a wider range of models for charged relativistic spheres. We illustrate this 
feature with an example involving a specific value for n. For example, suppose 
that n = 1 then h = (5 + a)a and we get 
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from (|2T|) where Oi and 02 are new arbitrary constants. From fl25l) and fl4aj) - 
fl4dp the solution to the Einstein-Maxwell system becomes 
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Note that the solution of the form fl26l) cannot be regained from Komathiraj 
and Maharaj [13] charged models except for the value of a 



(5+a) 



This indi- 
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cates that our model is the generahsation of Komathiraj and Maharaj charged 
models with more general behaviour in the gravitational and electromagnetic 
fields. 



6 Physical analysis 



In this section, we briefiy consider the physical features of the models generated 
in this paper. For the pressure to vanish at the boundary r = Rin the solution 
(I2U1) we require p{R) = which gives the condition 



4(1 + aCR 



2\2 



aa{h — a)CR^ 



(27) 



where y is given by (IT7|) - (ITU]) . This will constrain the values of a, h and a. The 
solution of the Einstein-Maxwell system for r > R is given by the Reissner- 
Nordstrom metric as 



2m q 



(28) 

where m and q are the total mass and the charge of the star. To match the 
potentials in (12(1 to (125]) generates the relationships between the constants 
Ai, y42, a, h and R as follows 



(l - ^ + =A'[Am{R) + A,y^{R)]\ (29a) 
2m g^r' 1 + bCR^ 

The matching conditions (12 7p and (1291) place restrictions on the metric coef- 
ficients; however there are sufficient free parameters to satisfy the necessary 
conditions that arise for the model under study. Since these conditions are 
satisfied by the constants in the solution a relativistic star of radius R is 
realisable. 

From (I20al) and (I20b[) we easily observe that the gravitational potentials e^'^ 
and e^'^ are continuous and well behaved for wide range of the parameters 
a and h. From (I20cp . the variable x can be expressed solely in terms of the 
energy density p as 
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Fig. 1. Energy density. 
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b'^C - 2bp - hc?C - 2abC ± ^(a-6)Cx 



X = 



2b{3a^C + bp) 



^[-27a?bC + a?C{2h + 66a) - h^ibC + 8p) + ab{?>bC + 8p + 2bap)] 

Hence, from ( ]20dl) the isotropic pressure p can be written as a function of en- 
ergy density p only. Therefore the solutions generated in this paper satisfy the 
barotropic equation of state p = p{p)- Many of the solutions found previously 
do not satisfy this desirable feature. We illustrate the graphical behaviour of 
matter variables in the stellar interior for the particular solution fl26l) . We 
assume that ai = —4.897, 02 = C = 1 and a = a = 1/4 for simplicity, and 
we consider the interval < r < 1. To generate the plots for p,p,E'^,dp/dp 
and p vs p, we utilised the software package Mathematica. The behaviour of 
the energy density is plotted in Fig. 1. It is positive and monotonically de- 
creasing towards the boundary of the stellar object. In Fig. 2, we have plotted 
the behaviour of matter pressure p, which is regular, monotonically decreasing 
and becomes zero at the vacuum boundary of the stellar object. In Fig. 3, we 
describe the behaviour of the electric field intensity. It is well behaved and 
a continuous function. In Fig. 4, we have plotted the speed of sound dp/ dp. 
We observe that < dp/ dp < 1 throughout the interior of the stellar ob- 
ject. Therefore the speed of the sound is less than the speed of the light and 
causality is maintained. In Fig. 5, we have plotted the pressure p verses the 
density p and we find that this approximates a linear function. This behaviour 
is to be expected as the gradients of p and p have similar profiles in the stel- 
lar interior. Thus we have demonstrated that the particular solution satisfies 
the requirements for a physically reasonable stellar interior in the context of 
general relativity. 



12 



p 




Fig. 2. Matter pressure. 

E- 

0.012 - ^ 



0.2 0.4 0.6 0.8 1.0 

Fig. 3. Electric field intensity. 
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